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Self-similar solutions for the dynamical condensation of a 
radiative gas layer 
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ABSTRACT 

A new self-similar solution describing the dynamical condensation of a radiative gas is 
investigated under a plane-parallel geometry. The dynamical condensation is caused 
by thermal instability. The solution is applicable to generic flow with a net cooling 
rate per unit volume and time oc p^T", where p, T and a are density, temperature 
and a free parameter, respectively. Given a, a family of self-similar solutions with 
one parameter ry is found in which the central density and pressure evolve as follows: 
p{x = 0,t) oc {tc - t)-''/(2-") and P{x = 0,t) (x [t^ - t)(i-'7)/(i-"), where <c is an 
epoch when the central density becomes infinite. For ~ 0, the solution describes the 
isochoric mode, whereas for ?7 1, the solution describes the isobaric mode. The self- 
similar solutions exist in the range between the two limits; that is, for < < 1. No 
self-similar solution is found for a > 1 . We compare the obtained self-similar solutions 
with the results of one-dimensional hydrodynamical simulations. In a converging flow, 
the results of the numerical simulations agree well with the self-similar solutions in 
the high-density limit. Our self-similar solutions are applicable to the formation of 
interstellar clouds (HI cloud and molecular cloud) by thermal instability. 
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1 INTRODUCTION 

Thermal instability (TI) is an important physical process 
in astrophy sical environments that are subject to radia- 
tive cooling. iFall fc Reed (|l985l l investigated the effect of TI 
caused by Bremsstrahlung on structure formation at galac- 
tic and subgalactic scales. TI also plays an important role 
in the interstellar medium (ISM). The ISM consists of a 
diffuse high-temperature phase (warm neutral medium, or 
WNM) and a clumpy low-temperature phase (cold neutral 
medium, or CNM), in supersonic turbulence. It is known 
that the ISM is thermally unstable in the temperature range 
between these the t wo stable phases; that is, in the range 
T ~ 300 - 6000 K (|Dalgarno fc McCravl[l97a '). Koyama & 
Inutsuka (2000, 2002) suggested that the structure of the 
ISM is provided by the phase transition that is caused by 

TI that occurs in a shock-compressed region. 

The basic properties of TI were investigated by iFieldl 

who investigated the TI of a spatially uniform gas 
using linear analysis for the unperturbed state in thermal 
equilibrium. The net cooling function per unit volu me and 
time was assumed to be Aop^T" — Fno. iFieldl (|l965l ) showed 
that gas is isobarically unstable for a < 1 and is isochori- 
cally unstable for a < 0. In thermal non-equilibrium, how- 
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ever, where the heating and the cooling rates are not equal 
in the unperturbed state, the stability condition of TI is 
different from that in thermal equilibrium. This situ ation is, 
for ex ample, generated in a shock-compressed region. iBalbusI 
1 19861 ) derived a general criterion for the TI of gas in a non- 
equilibrium state. A flow dominated by cooling is isobari- 
cally unstable for ot < 2 and is also isochorically unstable 
for a < 1. Therefore, the phase with T ~ 300 - 6000K is 
both isobarically and isochorically unstable. 

The non-linear evolution of TI has been inves- 
tigated by many a uthor s (iKovama fc Inutsukal l2002l : 
Audit fc Hennebelld l2005l: iHennebelle fc AuditI l2007l : 



Heitsch et al.l l2006l : IVazguez-Semadeni et ahF 2003) using 
multi-dimensional simulations. However, the nonlinear 
behaviour has only been investigated analytically by a 
few authors, using some simplifications in an attempt 
to ga in a deeper insight into the nature of TI. iMeersorJ 
(| 19891 ) investigated analytic solutions under the isobaric 
approximation. However, this approximation is valid only 
in the limit of short or intermediate-scale. 

In this paper, a new semi-analytic model for the non- 
linear hydrodynamical evolution of TI for a thermally non- 
equilibrium gas is considered without the isobaric approx- 
imation. Once TI occurs, the density increases drastically 
until a thermally stable phase it reached. During this 
dynamical condensation, the gas is expected to lose its 
memory of the initial condition, and is expected to be- 
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have in a self-similar manner (IZerdovich fc Raizeil 1 19671 : 
iLandau fc Lifshit j 1 19591 ). In particularly, in star forma- 
tion, self-similar (S-S) solutions for a runaway collapse of 
a self- RravitatinK gas have been i nvestigate by many au- 
thors jPenstonlll969l: lLarsonl | l969l: IShulll977l: lHunterlil977l : 



IWhitworth fc Summers! Il985l : iBoilv fc Lvnden-Belll Il995h . 
In this paper, S-S solutions describing dynamical conden- 
sation by TI are investigated by assuming a simple cooling 
rate without self-gravity. A family of S-S solutions which 
have two asymptotic limits (the isobaric and the isochoric 
modes) is presented. 

In Section (2] we derive the S-S equations, and describe 
the mathematical characteristic of our S-S equations as well 
as the numerical methods. In Section O our S-S solutions 
are presented and their properties are discussed. In Section 
|4l the S-S solutions are compared with the results of one- 
dimensional numerical simulations. The astrophysical impli- 
cations of the S-S solutions and the effect of dissipation are 
also discussed. Our study is summarized in Section (5] 



2 FORMULATION 
2.1 Basic equations 

In this paper, we consider S-S solutions describing runaway 
condensation at a region where the cooling rate dominates 
the heating rate. The S-S solutions describe the evolution of 
the fluid after its memories of initial and boundary con- 
ditions have been lost during the runaway condensation. 
Therefore, the solutions are independent of how a cooling- 
dominated region forms. 

A net cooling rate pjC = A — F is considered, where 
A and F are the cooling and heating rates per unit volume 
and time, respectively. In low-density regions, the main cool- 
ing process is spontaneous emission by coUisionally excited 
gases. In this situation, the cooling rate is determined by 
the collision rate, which is proportional to p^. Therefore, we 
adopt the following formula as the net cooling rate: 



p£{p,Ca) = Aop^C^" erg cm ^s \ 



(1) 



where P, Ca, 7 and a are the pressure, sound speed, the 
ratio of specific heat, and a free parameter, respectively. The 
effects of viscosity and thermal conduction are neglected for 
simplicity. 

The basic equations for a radiative gas under a plane- 
parallel geometry are the continuity equation. 



the equation of motion, 

1 d , 2 . „ 
Dt 7p ox 



and the energy equation, 
2 D In Cs D In p 



7-1 Dt 



Dt 



A 2(ct-l) 



(2) 



(3) 



(4) 



where D/Dt = d/dt + vd/dx indicates Lagrangian time 
derivative. 



2.2 Derivation of self-similar equations 

The density increases infinitely because the net cooling rate 
is assumed to be oc p^c^" in equation ([T]). The epoch at which 
the density becomes infinite is defined by tc- We introduce 
a similarity variable ^ and assume that the each physical 
quantity is given by the following form: 

x = xo{t)^, v{x,t) =vo(t)V{0, c4x,t)=co{t)X{^), 

p{x,t) = poitMO and P{x,t) = Po(t)n(C). (5) 

By assuming power- law time dependence of xo{t) = at", 
where t, = 1—t/tc and n = 1/(1 — a;), we find the following 
relations: 

xo{t) = at", uo(t) = co(t) = ^C", 

tc 



-20+3 „2(l-(i) 2 

-C", and Po{t) - P° " 



(6) 



Ao ■ ' 7 

where f3 = ~ 2a) — 1. Because the dimensional scale is 
introduced only by Ao , the time dependence of the physical 
quantities cannot be fixed. Therefore, we leave a; as a free 
parameter. 

For convenience, instead oi lv, we can use a parameter 
rj which is given by 

(2 - a){(3 - 2a)Lo - 1} 



V = 



1 - (X) 



(7) 



Using 77, the time dependence of pa{t) and Po{t) is given by 

po(t) oc t^/f^-"' and Po(t)octl^-'''/(i-"'. (8) 

Substituting equations (HI) and ^ into the basic equations 
((2|, Q and Q, the similarity equations are obtained as 



,,dlnf7 dV 



/3, 



7 



-ujV, 



(9) 
(10) 



and 



2 dlnX dlnO. 



7-1 de 



dC 



-uj - P ~ ■ynX^'-"-^'' . (11) 



2.3 Homologous special solutions 

The similarity equations ([9]), (|10p and (|lip have homologous 
special solutions with a spatially uniform density and sound 
speed. Substituting f2(5) — Qo and X(^) = Xo into the 
similarity equations ([9|, (|10p and (|lip . one obtains 

dC 

dV 

and 
2uj 



^_J-P- l^oX',^"-'^ = 0. 
Equation (|12|) is integrated to give 



(12) 
(13) 

(14) 

(15) 
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where V{(, = 0) = is assumed. Substituting equation (llSp . 
equation (|13|l becomes 



0. 



(16) 



Therefore, homologous solutions exist for /3 = 0, = or 
a — 1. Equation (|14|) gives the relation between fio and Xq. 
For l3 — 0, the solution is given by Sl(^) — ^o, V{0 — and 
X{^) — Xo- The time dependences of density and sound 
speed are given by 



p{t) = const, and Cs(t) cx tl^^^ 



(17) 



respectively. This is the homologous isochoric mode. For 
Q J5 1, the solution is unphysical because the sound speed 
is negative or complex. Therefore, no homologous isochoric 
solutions exist for a ^ 1. This fact is related to the Bal- 
bus criterion for the isochoric instability (|Balbusll 19861 ) . The 
solution with uj = indicates isothermal condensation. For 
a = 1, no S-S solutions are found except for the homologous 
solutions as shown in Section [31 



2.4 Asymptotic behaviour 

In this subsection, the asymptotic behaviour of the S-S so- 
lutions at ^ — > cxD and ^ ^ is investigated. 



2.4-1 Solutions for ^ —> oo 

Assuming that n(^) is a decreasing function of ^ and that 
I V(5)| and A (5) are increasing functions, equations ((91, (|10p 
and (|lip become 

dlnn dV^ dlnA 

^— ^—~ujV, and^-^~a; (18) 

to the lowest order. Expanding to the next order, the asymp- 
totic solutions are given by 

and 

A(e)~A^r + 

where 



(1) 



1 



1 - UJ 



LoVl + ^ {13 + 2cj) 



and 

(1) _ (7 + l)tjVooAoo + 7(7 - i)n^x, 
2(1 -w) 



2{q-1) 



(19) 
(20) 

(21) 

(22) 
(23) 

(24) 



2.4.2 Solutions for ^ ^ 

We expand S1(C), V{(,) and A(^) in the following forms: 



and 

XiO^Xoo+xg^e- 



(25) 



Substituting equations (|25p into ((9]), (jlOp and (jlll) . we obtain 
the following relations for coefficients; 



1/ « A 1/(1) "it^+'^M'o 

Voo = P and (,^,)^^^ ■ 



(26) 



The relation between Qqq and Aoo is given by 

7(7 - mooX^ot'''' = - /3(7 - 1). (27) 



The relation between fipg' and Aqq^ is given by 



(1) Y-(l) 
(lu 2 00 _ Q 

f^oo Aoo 



(28) 



and u is given by 
_ (2 - Q;)eoo 



where eoo = (7 - l)f^ooAo'<° '\ (29) 



2.5 Critical point 

Equations pop and (jlip can be rewritten as 



dlnO 
where 



72 

/i' 



/i = (v^ + 



dV 



/3 , dlnA 

— and 

h 



h 



X' 



9 + f, 



h 
/i' 



73 = -X^g + I3h -{V + 0.f, 



h = -- 



x''-j{v + 0' 



(7-l){f^ooAo^("-^) 



g(n,X) 
and 

f{^,v)=p{v + o~'^v. 



■ nA^'' 



(30) 

(31) 
(32) 
(33) 

(34) 
(35) 



A singular point exists where /i = in each of the equations 
(|30|l . This is the critical point (^ = ^s). In the S-S solutions 
which is smooth at this point, the numerators (72, I3 and 
74) must vanish at ^ = 0. Among four conditions Ii — I2 = 
13 = 14 = 0, there are two independent conditions which are 
given by 

A^ 



and 



v + c 



' + ,/ = 0. 



(36) 



2. 5. 1 Topological property of the critical point 

In this subsection, the topological properties of the critical 
point are investigated. We introduce a new variable s which 
is defined by ds = lid^ (jWhitworth fc Summerj|l985h . The 
dimensionless quantities at infinitesimal distances from the 
critical point, Ss -I- 5s, are given by 



Q{ss + 5s) = Qs+5Q, 



(37) 



where Q = (^ In f2, 1^, In A) . The subscript " s" indicates the 
value at the critical point. Substituting equation (|37p into 
equations (|30|) and linearizing, one obtains 



ds 



97, 



(38) 



The eigenvalues A of Aij provide the topological property of 
the critical point. The eigenequation is given by 
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■.X's) 



0, 



(39) 
(40) 



where 

61 = e, {(7 - 1)(q - 1) + 1} + + /3 - 2, 
and 

62 = 27(a - l)el 

+ {-Pi ~ 4(i)a + 2u) ~ 27a + 2q - 3 - 2/3 + a;7 + 7)es 

/? - 2a; + 2cj/3 



+ /3 + 2i^^ 



7 



(41) 



where 



(7 - i)ax, 



2{q-1) 



The degenerate solution, 



A = 0, is attributed to the homologous special solutions (see 
Subsection 12. 3|) . Therefore, the topology is described by the 
two eigenvalues other than A = 0. 



a=0.3 




a=0.5 




0.8 0.9 1 



a=0.8 




a=1.2 



SP 



0.8 0.9 1 

55 



= h, i = 1, 2, 3 and 4, 



(42) 



2.6 Numerical methods 

In this subsection, our numerical method for solving the 
similarity equations pop is described. We adopt s as an in- 
tegrating variable. Using s, equations pO|) are rewritten as 

dQj_ 
ds 

which are not singular at the critical point. 

The critical point is specified by (^s, Vs) because there 
are four unknown quantities Qs and two conditions, namely 
equations p6[l . Moreover, the similarity equations (|42p are 
invariant under the transformation ^ —> ?n^, V — > mV, 
n m~^^°'-^^n and X mX, where m is an arbitrary 
constant. Therefore, the S-S solution for (m^s,m14) is the 
same as that for (^s, Ki). As one can take m in such a way 
that (^s, V's) satisfies the relation + = 1, the critical 
point is specified only by ^s. 

Given ^s, the velocity is obtained as 14 = — \/ 1 — ^| . 
From equations (1361) . fis and Xs are determined. As Xs = 



Vs + £,B is positive, the range of is l/\/2 < < 1. The 
similarity equations (|42[) are integrated from to both di- 
rections (^ — > and ^ 00) along the gradient derived from 
equation (|38[) using the fourth-order Runge-Kutta method. 
The position of the critical point is determined itera- 
tively by the bisection method until the solution satisfies 
the asymptotic behaviours (see Subsection 12. 4p . 



3 RESULTS 

3.1 Physically possible range of parameters (17, a) 

In this section, we constrain the parameter space (77, a) by 
physical properties of the fiow. 

3.1.1 Range of rj 

Using equation ([8]), the time dependences of the central den- 
sity and pressure are given by 

poo(t) oc C"/'"-"' and Poo(t)cxil''''''''''"', (43) 

where the subscript "00" indicates the value at the cen- 
tre. During S-S condensation by cooling, the central density 
(pressure) must increases (decreases) monotonically with 
time. Therefore, the following two conditions are obtained: 



Figure 1. Topological properties of the critical point for a =0.3, 
0.5, 0.8, 0.9, 1.1 and 1.2. The ordinates and abscissas are 77 and 
^s, respectively. The letters 'N', 'SD' and 'SP' indicate a node, 
a saddle and a spiral, respectively. The open circles indicate S-S 
solutions. 



- 7]/(2 - a) < ^ ?7 > 
and 



(44) 



(l-77)/(l-a) >0 ^< ' , r 1 (45) 

^ "' ^ ' 1 r; > 1 for a > 1 . ^ ' 

From equations (|44|) and (|45p . the range of 77 is given by 

0<,<1 for a<l 

77 > 1 for a > 1 . ^ ' 

For 77 ~ 0, the time dependence of poo and Poo is given by 

Poo (t) ~ const, and Poo(t) oc ti''*^""^ (47) 

This time evolution indicates the isochoric mode. For 77 ~ 1, 
the time dependence of poo and Poo is given by 



poo(t) oc "^'''^ °' and Poo(t) — const. 



(48) 



corresponding to the isobaric mode. The critical value of 
r\ = '7oq is given by the condition that the increasing rate 
of poo(i) is equal to the rate of decrease of Poo(t), which is 
given by 



r7oq = (2-a)/(3-2a). 



(49) 



For 77 = 77eq, the time dependences of poo and Poo are given 
by 



poo 



(t) oc and Poo(t)octy(^-^'^). 



(50) 



Therefore, the S-S solutions for < 77 < 77oq and 7ycq < 
77 < 1 are close to the isochoric and the isobaric modes, 
respectively. 



3.1.2 The range of a 

In the similarity equations (|42[) . the critical point is specified 
by for given a and 77. The topology of the critical point is 
given by the two eigenvalues of the eigenequation (|39p other 
than A = 0. When the two eigenvalues are real and have 
the same sign, the critical point is a node. When the two 
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Figure 2. The distribution of dimensionless quantities {Q,—V,Tl,X) for (a)r;=0.98, (b)0.76 and (c)0.12. The open circles indicate the 
critical points. Panels (d) and (e) indicate the dependence of Q/Qoo and n/IIoo, respectively, on rj. 




Figure 3. Time evolution of dimensional quantities (p, P, —v). The solid, long-dashed, short-dashed and dotted lines are corresponding 
to t,=l, 0.1, 0.01 and 0.001, respectively. The columns correspond to »?=0.98, 0.76 and 0.12 from left to right, respectively. The arrows 
indicate the time evolution of flow. 



eigenvalues are real and have the opposite signs, the critical 
point is a saddle. When the two eigenvalues are complex, 
the critical point is a spiral. Fig. [T] shows the topology of 
the critical point in the parameter space (^s, 77) for a =0.3, 
0.5, 0.8, 0.9, 1.1 and 1.2. In Fig.[Tl the letters 'N', 'SD' and 
'SP' indicate a node, a saddle and a spiral, respectively. The 
open circles indicate the numerically obtained S-S solutions 
presented in Subsection 13.21 Fig. [T] shows that the topology 
drastically changes at a = 1. For a < 1, because the critical 
point is a node or a saddle, the S-S solutions are expected to 
exist. For a ~ 1.1, a non- negligible fraction of the parameter 
space is covered by a spiral. Almost all parameter region is 
a spiral for a > 1.2. We searched the S-S solutions that 
connect < ^ < 00 and found them only in N and SD for 
Q < 1. 



3.2 Self-similar solutions 

In this subsection, we describe new S-S solutions and inves- 
tigate their properties. As shown in Subsection 13.1.11 our 
S-S solutions exist for < < 1 and a < 1. 

FigsEJa), (b) and (c) show the obtained S-S solutions 
for a = 0.5 with 77=0.98, 0.76 and 0.12, respectively. The or- 
dinates denote the dimensionless quantities {^l, —V, 11, X). 
The abscissas denote the dimensionless coordinate ^. The 
S-S solutions for 77=0.98 and 0.12 approximately represent 
the isobaric and the isochoric modes, respectively. The so- 
lution for 77 = 0.76 ~ 77cq corresponds to the intermediate 
mode shown in Subsection 13.1.11 The open circles indicate 
the critical points. In Figs[2ja), (b) and (c), one can see that 
the S-S solutions depend strongly on -q. Let us focus on the 
dependence on rj of density and pressure, shown in Figs[2][d) 
and (e), respectively. In these figures, it can be seen that, for 
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n 


-r?/(2 - a) 


(l-r,)/(l-a) 


€s 


Hoc) 


^00 








0.059 


1.923 


-0.039 


1.882 


0.9999 


0.873 


0.977 


0.882 


-0.0234 


0.0408 


0.272 


1.818 


-0.181 


1.455 


0.9975 


0.768 


0.883 


0.782 


-0.124 


0.201 


0.500 


1.667 


-0.333 


1.000 


0.9907 


0.643 


0.766 


0.612 


-0.244 


0.391 


0.761 


1.493 


-0.507 


0.615 


0.9792 


0.486 


0.609 


0.343 


-0.333 


0.613 


0.945 


1.370 


-0.631 


0.110 


0.9853 


0.360 


0.472 


0.123 


-0.223 


0.805 


0.995 


1.337 


-0.663 


0.011 


0.9986 


0.237 


0.315 


0.0260 


-0.0576 


0.952 



Table 1. Relevant parameters for obtained self-similar solutions (q = 0.5) are summarized. 



T] = 0.98 (the isobaric mode), the density increases sharply 
towards ^ = 0, whereas the pressure is almost spatially con- 
stant. In contrast, for rj — 0.12 (the isochoric mode), the 
density is almost spatially constant whereas the pressure 
decreases sharply towards ^ = 0. 

Next, the S-S solutions are shown using the dimensional 
physical quantities for various rj. Fig. [3] shows the time evo- 
lution of p (first row), P (second row) and —v (third row). 
The columns in Fig. [3] corresponds to ry=0.98, 0.76 and 0.12 
from left to right. The solid, long-dashed, short-dashed and 
dotted lines correspond to t*=l, 0.1, 0.01 and 0.001, respec- 
tively. From Fig. [S] one can clearly be seen that 77 specifies 
how the gas condenses. In the case with rj = 0.98 (the iso- 
baric mode), runaway condensation occurs and the density 
increases in the central region, whereas the pressure is spa- 
tially and temporally constant. In the case with 77 — 0.12 
(the isochoric mode), the density is spatially and tempo- 
rally constant and Poo decreases considerably. In the case 
with ?7=0.76, the condensation occurs in an intermediate 
manner between the isobaric and the isochoric modes. 

Relevant parameters, {rj, n, —rj/{2~a), (1 — 77)/(l — q), 
^s, Qoo, Xqo, Qao, Vac, ^oo), for the S-S solutions obtained 
for a = 0.5 are summarized in Table [T] 

From the above results, it is confirmed that the S-S so- 
lutions are specified by rj and contain two asymptotic limits 
(the isobaric and the isochoric modes). Moreover, the pa- 
rameter rj is related to the ratio of the sound-crossing to the 
cooling time-scales in the central region. Both time-scales 
are defined by 



t 



and 



i(^) ^ccn 



COO (t) nXoo 



Poo(t) 



(7- l)Aopoo(t)2coo(f)2 



^-1 y2{l-Q) 

"00 ^00 

7(7-- l)n 



tct* 



(51) 



(52) 



respectively, where cqo is the sound speed at the centre, and 
the scale length Xccn{t) and ^con are defined by 



p{x,,n{t),t)/p{0,t) = Q(Cccn)/f^OO = 0.8. 



(53) 



From equations (|51|l and (|52[) . the ratio r of the above two 
time-scales is obtained as 



''CI 



^con 
cool 



■ ?con7(7 ~ I)fio0-fo0 ^• 



(54) 



Fig. |4] shows the dependence of r on rj for a=-0.5, 0.5 and 
0.8. The figure shows that r ^ for rj ^ 1 and r — > oo for 
?7 ^ independent of a. For ~ 1, because t <^ 1, the gas 
condenses in pressure equilibrium with its surroundings. For 
?7 ~ 0, because r ^ 1, the pressure decreases and the density 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Figure 4. The dependence of r (equation I54| l on rj. The open 
circles, filled circles and open boxes indicate the S-S solutions for 
a=-0.5, 0.5 and 0.8, respectively. 



increases only slightly. These facts are consistent with the 
discussion in Subsection 13.1.11 and with Fig. [3l 

Finally, we consider the range of a for which the S-S 
solutions exist. When a = 1.0, no S-S solutions are found 
except for the homologous special solutions (see Subsection 
12. 3[) . Moreover, as mentioned in Subsection 12.31 there are 
no homologous isochoric solutions for a > 1. Therefore, for 
a > 1, the condensation of the fluid is not expected to be 
self-similar. iFouxon et all (|2007l ) found a family of exact so- 
lutions for Q = 1.5 > 1. Their solutions are not self-similar 



3/2 



lowing different scaling laws: igoSnd oc 
In their solution, during condensation, as tg, 



becomes 



expected to be constant and the isobaric approximation is 
locally valid. 



4 DISCUSSION 

4.1 Comparison with the results of 

time-dependent numerical hydrodynamics 

We compare the S-S solutions with results of time- 
dependent numerical hydrodynamics using one-dimensiona l 
Lagrangian 2nd-order Godunov method (jvan Leedll997l ). 
A converging flow of an initially uniform and thermal- 
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Figure 5. The time evolution of poo (thick Unes) and Poo (dashed lines) as a function of 1 — t/tcon for (a)L=0.01, (b)0.02, (c)0.2 and 
(d)0.5. The thin solid lines indicate the corresponding S-S solutions. 
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Figure 6. The distribution of the rescaled density p(l — t/tcon)'"''^ "'^ and pressure P(l — t/tcon) °') as a function of scaled 

coordinate x{l ~ t/tcon)~" ■ The L and rj in each panel are the same as that in Fi§. ^\ The lines are labelled according to 1 — t/tcon- 
The rescaled density and pressure at the centre are set to 10 and 0.1, respectively. The thin solid lines indicate the corresponding S-S 
solutions. 



equilibrium gas {p{x) = P{x) = 1), is considered. An ini- 
tial velocity field is given by v{x) — — 2 tanh(a::/L), where 
L is a parameter. The net cooling rate is assumed to be 
pC{p,T) = {p^/'^P^''^ - p)/{-i-l), where a=0.5. The inflow 
creates two shock waves, which propagate outward. Run- 
away condensation of the perturbation owing to TI occurs 
in the post-shock region because the cooling rate dominates 
the heating rate. The central density continues to grow and 
ultimately becomes infinite. During the runaway condensa- 
tion, the flow in the central region is expected to lose the 
memory of the initial and the boundary conditions, and to 
converge to one of the S-S solutions. The scale of perturba- 



tion in the post-shock region can be controlled by the value 
of L. The larger L is, the larger the scale of perturbation be- 
comes. We perform numerical simulations for the case with 
L=0.005, 0.01, 0.02, 0.2 and 0.5. 

Fig. [5] shows the time evolution of pcon(t) and Pcsnit) 
as a function of 1 — t/tcon, where pccn and Pccn are the cen- 
tral density and pressure, and icon is an epoch when pcen 
becomes infinity estimated from the numerical results. The 
panels represent (a)L=0.01, (b)0.02, (c)0.2 and (d)0.5, re- 
spectively. Using equation (|43|) . the time evolution of pccn(t) 
and Pcen (i) of each simulation provides the corresponding 77. 
As a result, the corresponding values of 77 are found to be 
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(a) 0.99, (b)0.98, (c)0.76 and (d)0.50. The thin soUd Unes 
in Fig. [5] indicate the corresponding S-S solutions. From the 
value of T] and Fig.O cases (a) and (b) correspond to the iso- 
baric mode, case (c) corresponds to the intermediate mode 
iv ~ rjcq) and case (d) corresponds to the isochoric mode. 
There is a difference in the convergence speed from the cor- 
responding S-S solution. Fig. [5] shows that the convergence 
speed of (a) and (b) is faster than that of (c) and (d). This is 
because the scale-length of the condensed region of (a) and 

(b) is smaller and the sound-crossing time-scale is shorter. 
Fig. [6] shows the time evolution of the rescaled density 

p(l - t/tcon)'"''^^"' and pressure P(l - t/tcon)"'^""'/'^""' 
as a function of the rescaled coordinate x{l — t/tcon)~" ■ The 
parameters L and 77 in Fig. [6] are the same as those in Fig.[5l 
The thin solid lines in Fig. [S] indicate the corresponding S-S 
solutions. In all the cases (a-d), the central regions are well 
approximated by the corresponding S-S solutions. From Fig. 
[S] it can be seen that the region that is well approximated 
by the S-S solution spreads with time in each panel. 

From the above discussion, the non-linear evolution of 
perturbation is well approximated by one of the S-S solu- 
tion in high density limit. Which S-S solution (0 < 7; < 1) is 
more likely to be realised in actual environments where var- 
ious scale perturbations exist? To answer this question, the 
dependence of icon on L is investigated. Fig. [7] shows the de- 
pendence of 1/fcon on 1/L. Here, 1/L and 1/tcon correspond 
to the wavenumber and the growth rate of the perturba- 
tion, respectively. Therefore, Fig. [7] can be interpreted as 
an approximate dispersion relation for the non-linear evo- 
lution of TI. For larger 1/L, 1/fcon is larger and 1/tcon 
appears to approach a certain asymptotic value. A similar 
behaviour is also found in the dispersion relation of iFieldl 
(| 19651 ) without heat conduction. In Fig. [71 the condensation 
for 1/L — 5.0 corresponds to the intermediate mode. Fur- 
thermore, although the isobaric condensation grows faster, 
there is little difference between 1/tcon for 1/L=200 (the 
isobaric mode) and 2 (the isochoric mode). Therefore, both 
the isobaric and the isochoric condensations are expected to 
coexist in actual astrophysical environments. 



4.2 Astrophysical implication 

In astrophysical environment, there are roughly two ranges 
of temperature with a < 1 where our S-S solutions are ap- 
plicable. 

In the high-temperatu re case, T > 2 x lO^K 
ISutherland fc Dopital 1 1 9931 ). the dominant coolants are 
metal lines and Bremsstrahlung. In a shock with a veloc- 
ity greater than 

^^120kms-(^^)''^ (55) 

the post-shock region is thermally unstable (a < 1). For ex- 
ample, a typical supernova satisfies this condition, and the 
S-S condensation is expected to occur. In reality, the den- 
sity does not become infinite because the cooling time-scale 
increases when the gas reach the thermally stable phase. In 
this case, the cooled layer will rebound and overstability may 
take place ([Chevalier fc Imamural Il982l : lYamada &: Nishil 
I2OOII : [Sutherland et al.ll2003l ). 

In the low-temperature case, 300< T <6000K, the dom- 
inant coolant is neutral carbon atom. iKovama fc Inutsukal 



1.7 
1.6 
1.5 
1.4 
1.3 
1.2 
1.1 



— 1 — I — I I I I 1 1 1 1 — I — I I I I 1 1 



— 1 — I — I I I I 




_i I ' I I I ' 



—I I ' ' ' I 



10 100 
1/L 



1000 



Figure 7. The dependence of the growth rate, 1/tcon, on the 
wavenumber, 1/L. 



(|2000l ') suggested that the structure of the ISM is caused by 
TI during the phase transition (WN M CNM). The phase 
transition is induced by a shock wave. lHennebelle fc Perauld 
(jl999f l considered the converging flow of a WNM under a 
plane-parallel geometry using one-dimensional hydrodynam- 
ical calculations, and investigated the condensation condi- 
tion of the Mach number and pressure in the pre-shock re- 
gion. In situations which satisfy their condition, S-S conden- 
sation is expected to occur. 



4.3 Effects of dissipation 

In the actual gas, the effect of viscosity and heat conduc- 
tion becomes important on small scales. The importance of 
dissipation can be evaluated by the ratio between the advec- 
tion and the dissipation terms, or Raynolds number, which 
is given by 



7^ = 



Cvpvd^T 
d4K{T)d^T) 



c 

K{Too) 



(56) 



where Cv, K{T) and U are the specific heat at constant vol- 
ume, the heat conduction coefficient and the typical veloc- 
ity, respectively. The Plandtl number is implicitly assumed 
to be of order unity. A typical length-scale is assumed to be 
Xccn{t). If the flow converges to one of the S-S solutions, 
the ratio, r, between the sound-crossing and the cooling 
timescales is constant. Therefore r is defined by 



^sound _ _ / 1 \ A „ „2a — 3 

— - = a::con7(7 " l)AopooCoo . 



Using T, Xccn is written as 



-1 S-2a 
Poo C()0 



7(7 - l)Ao 



(57) 



(58) 



Substituting this equation into equation (|56p . one obtains 

CvfvT 4-2Q-K 



7^ 



7(7 - i)Ao/«:o 



(59) 



where U = fvCooifv < 1) and K{Too) = KocSo- 

For the low-temperatur e case, we a dopt K — 2.5 x 
10^ VT ergs cm "-^ s~^ (|Parkeij|l953l') a nd the following 
cooling function (jKovama fc InutsukalbOO^ ): 
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A(r) ~ 1.0 X 10^° p^Vt erg cm"^ s"\ (60) 
From equation (|59|) . the Raynolds number is given by 

In the isobaric mode, r ~ 0.1 from Fig. U Because the 
Raynolds number is much larger than unity, the dynami- 
cal condensation of the post-shock region is expected to be 
well described by the S-S solution with r > 0.1. 

For the high-temperature case, the cooling rate of metal 
lines (lO^K < T < lO^K) is given by 

A(r) ~ 2.2 X IO'V'T""-® ergs cm"' s"\ (62) 
We a dopt K = 1.24 x 10-6^5/2 ^^^^ ^^-1 ^-1 g-i (|p^^^ 
I1953I ). Using this formula and eqution (|59[) . the Raynolds 
number is given by 

^-«Klr)((&)(T(?K)"- 

Because of 7?. 3> 1, in the high-temperature case, the S- 
S solution with r > 0.1 is expected to well describe the 
dynamical condensation. 

As the temperature decreases, TZ decreases and becomes 
less than unity at a certain epoch. Thereafter, our S-S so- 
lutions are invalid and the effect of dissipation becomes im- 
portant. 

5 SUMMARY 

We have provided S-S solutions describing the dynamical 
condensation of a radiative gas with A — F ~ p^T". The 
results of our investigation are summarized as follows: 

(i) Given a, a family of S-S solutions written as equation 
([8]) is found. The S-S solutions have the following two limits. 
The gas condenses almost isobarically for 77 ~ 1, and the gas 
condenses almost isochorically for 77 ~ 0. The S-S solutions 
exist for all values of 77 between the two limits (0 < < 1). 
No S-S solutions are found in a > 1. 

(ii) The derived S-S solutions are compared with the re- 
sults of numerical simulations for converging flows. In the 
case with any scale of perturbation in the post-shock re- 
gion, our S-S solutions approximate the results of numerical 
simulations well in the high density limit. In the post-shock 
region, smaller-scale perturbations grow faster and asymp- 
totically approach to the S-S solutions for ~ 1. However, 
the difference of the growth rate between the isobaric mode 
iVcq < 77 < 1) and the isochoric mode (0 < < 77cq) is 
small. Therefore, in actual astrophysical environments, the 
S-S condensations for various 77 are expected to occur simul- 
taneously. 

In the S-S solutions presented in this paper, a plane- 
parallel geometry is assumed. This situation is much simpler 
than the two-dimensional simulations that are currently be- 
ing carried out involving TI, so the results are of limited 
use, although the results may be relevant to code tests that 
are run on a one-dimensional simulation. The stability anal- 
ysis on the present S-S solution may provide some insight 
into the multi-dimensional turbulent velocity field. The ef- 
fect of magnetic field and self-gravity may also be important 
in ISM but they are neglected in this paper. These issues will 
be addressed in forthcoming papers. 
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